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Abstract: We extend the study of corrections to the eikonal approximation that was 
initiated in Ref. [1] to higher orders. These corrections associated with the finite width of 
the target are investigated and the gluon propagator in background field is calculated at 
next-to-next-to-eikonal accuracy. The result is then applied to the single inclusive gluon 
production cross section at central rapidities and the light-front helicity asymmetry, in pA 
collisions, in order to analyse these observables beyond the eikonal limit. The next-to- 
next-to-eikonal corrections to the unpolarized cross section are non-zero and provide the 
first corrections to the usual &_L~f ac t° r i ze d expression. In contrast, the eikonal and next-to- 
next-to-eikonal contributions to the helicity asymmetry vanish, while the next-to-eikonal 
ones are non-zero. 
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1 Introduction 


Particle production in high-energy collisions involving nuclei constitutes a hot topic of 
research in Quantum Chromodynamics (QCD) due to both its theoretical and practical 
importance. Depending on the size of the relevant scale associated to the production 
process, extensions of collinear factorization (see [2] and references therein) for a truly 
hard scale, or of fcj_-factorization [3, 4] for a semi-hard scale, are demanded. 

The existing calculations for single inclusive particle production in pA collisions - a 
situation usually referred to as dilute-dense scattering - result in a /cj_-factorized expression 
[5-10] at the lowest order in the coupling constant if the particle is produced at central 
rapidity. By contrast, the hybrid factorization [11] is relevant if the particle is produced 
at forward rapidity 1 . All these calculations rely on: (a) the replacement of the target 
by a large background field; and (b) the use of the eikonal approximation, in which the 
constituents of the projectile just experience color rotation upon scattering with the target 
through picking a Wilson line at a given transverse point but integrated along the light- 
cone direction of propagation of the projectile, for which the target appears as infinitely 
Lorentz contracted. The Color Glass Condensate (CGC, see [19] and references therein) 
employs both assumptions for resumming high-energy logarithms, but neglects corrections 
suppressed by the energy of the collision, a situation parallel to neglecting the inverse 
powers of the hard scale in collinear factorization. 

In the previous work [1] we developed a systematic method to compute higher order 
corrections to the eikonal approximation in the CGC (see related work in the context of 
soft gluon resummation for hard scattering amplitudes in [20, 21]). Such next-to-eikonal 
contributions, due to the non-zero length of the target or the finiteness of the energy 
of the projectile, can be understood as the subleading effects with respect to either the 
infinite Lorentz contraction of the target or the infinite Lorentz dilation of the projectile. 
The method involves new operators that contain transverse gradients of the background 
field in the target, which are thus sensitive to its transverse structure. We applied such 
formalism to single-inclusive particle production at central rapidities, recovering the usual 
k ±-factorized formula in pA collisions and obtaining the first next-to-eikonal correction. 
For collisions in which no privileged direction exist, the next-to-eikonal correction was 
found to vanish identically, thus extending the validity of the usual fcj_-factorized formula. 
On the contrary, for collisions in which there is some privileged direction like in the study 
of spin asymmetries for a polarized target, it was the next-to-eikonal correction the first 
non-vanishing contribution. This resembles the situation in collinear factorisation where 
the first power-suppressed contributions become dominant for semi-inclusive production or 
polarized collisions, see e.g. [22-28] and references therein. 

In the present paper we extend the formalism of [1] to next-to-next-to-eikonal accuracy. 
In Section 2 we compute the corresponding retarded gluon propagator that will be used 

1 There exists a large activity on computing and understanding the corresponding next-to-leading cor¬ 
rections in the coupling constant [12-18]. 
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subsequently to compute particle spectra. In Section 3 we calculate both unpolarized and 
polarized single inclusive gluon production in unpolarized pA collisions. Finally we present 
our conclusions, that we anticipate here: the computed next-to-next-to-eikonal corrections 
for the gluon helicity asymmetry vanish identically, while those for unpolarized production 
are non-zero and provide the first corrections to the fcj_-factorized expression [5-10]. These 
results extend the conclusions extracted in the previous work [1] to higher orders. New 
operators appear, whose meaning is discussed and whose rapidity evolution is left for future 
studies 2 . 

2 Retarded gluon propagator in background field at next-to-next-to- 
eikonal accuracy 

In the calculation of high-energy dilute-dense scattering processes, for example within the 
CGC formalism, the most important building blocks are the gluon and quark propagators 
in a strong classical background gluon field representing the dense target. In this paper 
we restrict ourselves to the gluon propagator. We choose a frame in which the target is a 
highly boosted left-mover, and pick the light-cone gauge, A+ = 0. 

In the limit of large Lorentz boost, one can in principle perform the following approx¬ 
imations to the classical background field Aa(x): 


4£(®) - ^ 4 a (x), 

(2.1) 

4£(*) - 40(z + ,x), 

(2.2) 

A%{x) oc d(x + ). 

(2.3) 


These approximations are together equivalent, at the level of the background field, to the 
eikonal approximation [33, 34] usually formulated at the level of the scattering amplitude. 
The second and third approximations correspond respectively to infinite time dilation and 
infinite length contraction. Corrections to each of the three approximations (2. 1-2.3) are 
suppressed by inverse powers of the Lorentz boost of the target, which correspond to inverse 
powers of the total energy at the level of a high-energy scattering amplitude for a dilute 
projectile off that background field. 

In the present paper, as in Ref. [1], we investigate corrections to the eikonal approx¬ 
imation related only to the violation of the approximation (2.3), i.e. associated with the 
finite length of the target. Although the three approximations (2. 1-2.3) have a priori the 
same validity range, it can be justified to relax only the approximation (2.3) when the 
target is a large nucleus, since the A 1 / 3 enhancement of the length can compensate the 
suppression by the boost factor. 

Hence, in the rest of the present section, we will study the retarded gluon propagator 
G 1 ]^ ( x , y) a b in the presence of a classical gluon background field 

4£(®) = A“(x + ,x), (2.4) 

with an arbitrary x + and x dependence. 

2 Note that the rapidity evolution of these decorated dipole operators may be related to the resummations 

recently discussed in the context of jet quenching in Refs. [29-32]. 
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2.1 From gluon to scalar propagator in the background field 

In the light-cone gauge, the gluon propagator obeys the constraints 

Gn'{x,y)ab = G^~(x,y) ab = 0. (2.5) 


Due to the x independence of the background field (2.4), it is useful to take the Fourier 
transform from x~ to k + , and define 

G^,«u = f ( 2 . 6 ) 

using the notation x = (x + , x). 

Linearizing the classical Yang-Mills equations around the background field M^(x + ,x), 
one can find the Green’s equations obeyed by the background propagator (x, y) a b- 
Using the Fourier representation (2.6), the solution of these Green’s equations can be given 
in terms of a retarded scalar propagator in background field, as 


Gt(x;y) ab = 5 ij gfi-fay), 

£*:+ (x; y) ab = $k+ fe v )» 

^Ax-,yJ ab = ^ e d r^Ax-,yJ, 

$k+ (s y) ab = ^ 5 yi Qt + fe y) + ^\x-y ), 


(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


where the scalar propagator (x; y ) is the retarded solution of the Green’s equation 


safe 


id. r + + 


d 2 

2(fc++ie) 


+ g( y A~(x)-Tj b g b k c + (x-,y)=iS ac 6W(x-y). (2.11) 


We refer the reader to Ref. [1] for more details on the procedure just described. 

The interpretation of the solution (2.7-2.10) for the gluon propagator in background 
field is the following. Due to the approximations (2.1) and (2.2), interactions with the 
background field M~(:r + ,x) cannot change the polarization state of a gluon. Hence, the 
two physical polarizations of the gluon propagate independently through the medium, 
each of them according to the scalar equation (2.11). Then, the additional term with 
5 ( 3 )(x —y) in eq. (2.10) accounts for the non-propagating polarization associated with 
purely Coulombian interaction. 

The scalar Green’s equation (2.11) has a form of a Schrodinger equation in 2 + 1 
dimensions with a matrix potential. Its solution can thus be written as a path integral [35- 
39]. For our purposes, we need the discretized form of the path integral, which reads 





"TV—1 

| £o,fc+(4+l> Z n+U^n> Z a) 

_n= 0 


U ab (x + ,y + ;{z n }) , 

( 2 . 12 ) 
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with the boundary conditions zo = y and zjy = x, and 

4 = V + + ^ + -y + ). (2.13) 

The free scalar propagator appearing in the path integral expression, given in eq. (2.12), 
reads 

S 0M (x,y) = 9(x + -s, + ) ( 2jr ~f_ + g+) ) e*P ( 2 (x‘-y+) (X - y);i ) ' <2 ' 14) 


The discretized Wilson line is defined as 


U ab (x + ,y+,{ z n }) 



N-l 
II exp 

n=0 


ig 





(2.15) 


with V+ denoting path ordering along the .x + direction. 

In the expression of observables related to dilute-dense scattering at high energy, 
the background propagator appears typically through its transverse Fourier transform 
f d 2 x e _ * k x Q%+{x,y). In order to study corrections with respect to the eikonal approx¬ 
imation (or more precisely, with respect to the shockwave approximation (2.3)) for the 
background gluon propagator, we start from that Fourier representation of the background 
scalar propagator, and define the corresponding medium modification factor 7 Z c /. b (x + , y + ; y) 
as 


'2^ r ik-x nab , 


crx e 


gf + (x,y) = nt b (x+,y+- y) / d 2 x e 




^0 ,k+ 


(*, v) 


= ni b {x + ,y + ;y) 6{x + -y + )e 


-zk-y 


— — k 2 

> 2 fc+ 


(2.16) 


In order to calculate the expansion of TZ^. b (x + , y + ; y) beyond the eikonal approximation, 
it is more convenient in practice to consider the large k + limit than the limit of large boost 
of the target. The two limits are of course equivalent, thanks to the invariance of the whole 
scattering process under longitudinal boosts. 

There are two ways to define the large k + limit: either it can be taken at fixed k, or at 
fixed k/£: + . The former case is relevant for high-energy scattering processes in the Regge 
limit, whereas the latter is relevant for hard scattering processes in the Bjorken limit. The 
variable k /k + is related to the deflection angle for the gluon. One considers small angle 
scattering in the first case and finite angle scattering in the second, in the large k + limit. 

Both for convenience, and in order to obtain results applicable in each case, we perform 
the large k + expansion of 1Z^ b (x + , y + ; y) in two steps. First, we calculate the k + expansion 
at fixed k /k + in the section 2.2, and then we re-expand the result for small k /k + in the 
section 2.3. 


2.2 Semi-classical expansion for the background propagator 

In the large k + limit at fixed k/fc + , the path integral for the gluon propagator becomes 
increasingly dominated by the classical trajectory 

z(z + ) = y + (z + -y + )^~. (2.17) 
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Our method to calculate the expansion of lZ ( £ > {x+,y+] y) in this limit is the following. 
First, we take the Fourier transform x i— > k of the discretized expression (2.12) for the 
propagator and change variables in order to write the path integral as an integral over 
deviations with respect to the classical path (2.17), making evident the factorization of the 
free propagator contribution as in eq. (2.16). Second, we take the continuum limit. Third, 
we Taylor-expand insertions of the background field around the classical path, and perform 
explicitly the path integration order by order, in order to get the eikonal expansion at finite 
angle. 


2.2.1 Discretized form of the scalar medium modification factor 

Taking the Fourier transform x i —> k of eq. (2.12) and relabeling x —> z tv, one gets 


N 


d 2 x e- ?:k x g$(x,y) = Jin^ 


x e 


—ik-Zjv 


N -1 

n 

_n=0 


exp 


ID 2 

\n=l 

ik+N 


-ik+N 


N 


2 (x+ -y+) 


2ir(x+ - y+) 

(z n +i - z n ) 2 ^ U ab {x+,y+\{ z„}) 


(2.18) 


where zq = y. The next step is to perform the change of variables z n H > u n = z n — z n , 
where z n is the discretized version of the classical trajectory, i.e. 

k 


Zn = y + (z+ - y+) 


k+ 


Noting that 

(Zn+l Z n ) — 


k 

] k+ 


("n+1 “ Z t )tN + U «+l “ U « 


^(x + ~y + )j^:+ u n+l - U„ 


(2.19) 


( 2 . 20 ) 


= (u n+ i - u n ) 2 + — (x + - y + ) — ■ (u n+ i - u n ) + 


(x+ - y + ) 2 / k 


one finds that 


a —ik-zjv 


N -1 

n 

_n =0 


exp 


N 


ik+N 


2(x+ - y+) 


k+ 


(Zn+l z r 


N 2 \k+ 


k . ( x ' — y 1 )k •i 

■y - i ~—^- ,ku v 


N -1 

n ex p 

_n=0 


= P -* k y 


ik+N 


2(x+ - y+) 

N-l 


(Un* 1 U„ 


Jk- (uat-uo )+i ( * + 2fc +~' )k2 


(3!+-l/+)k- 


n 

_n=0 


exp 


ik+N 


2(x+ -y+) 


(u n +i 


( 2 . 21 ) 


because uo = 0. From the equations (2.16), (2.18) and (2.21), one obtains the discretized 
path-integral expression for the medium-modification factor IZ^ (x+, y+] y) as 


N 


N-l 


lZf{x+,y+-,y) = \im^ j ( \\d 2 u n j <j g oM (z+ +1 ,u n+1 ; z+ ,u r 


\n=l 


n =0 


xU ab (x + , y+ ; {z n + u n }) 


( 2 . 22 ) 
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2.2.2 Continuum limit 

Let us remind the Lie product formula which is valid for any square matrices A and B: 


e A+B = lim feT-efV. (2.23) 

T-J. + 00 \ ) 

For small B, eT = 1 + 0(B). Thus, for small B, the Lie product formula can be written 
as 


„ A+B 


= i lim + |e‘ — lj e '' (2.24) 

3 =i 

i -1 l-j 

' — lj e * (e* — lje * + 0(B 3 

3 =1 i=l 


l —^-(-OO 


(i-i)A 


( l-j-i)A 


In the large l limit, each sum can be written as an integral and one obtains 

e A+B = e A + f ds e sA B e (l-a)A + f ds f du e *A B e uA B g (l -s-u)A + ()(&). 

Jo Jo Jo 

(2.25) 


Introducing the notations 


a n = (x + -y + ) igT ■ A (z+; z n ) 
b n ( u«) = (x + -y + ) igT ■ [A~(z+; z n + u n ) - A~(z+] z n )\ , 


(2.26) 

(2.27) 


and using eq. (2.25), one can express each elementary gauge link (see eq. (2.15)) of the 
discretized Wilson line appearing in eq. (2.22) as 


exp 


(x+-y+) . 


N 


igT-A (z+;z„ + u n ) 


= e n 


»l 


= e n + / ds e s n 

Jo 


Bn. b n (u n ) 
n —-—- t>y ± *) n 


N 


(2.28) 


+ 


if j s— ^n(Wi) u— ^n(Uti) (1—s—u) — 

ds / due s N —~—— e * —V— e u u ) N 


N 


N 


+ O 


bn (Un) 

N 


This expansion can be inserted at each discrete step in eq. (2.22). 

In order to understand what happens in the continuum limit N —> +oo, let us first 
consider the case in which only the leading term e is kept at each step in eq. (2.22), 
which we call the contribution 0 to TZk(x + , y + ] y). One has 


Kk(x + ,y + ] y) 


= lim 

TV—>oo 


' N \ TV-1 ( \ 

d U n j T + < @0,k+ ( Z n+1 ’ u «+l) z n " u n) e N f 

,n=l J n=0 l J 


= lim U(x + ,y + ,{z n }) 
TV—>oo 

= U ( x + ,y + ; [z(z + )]) , 


(2.29) 
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using the identity 


J d 2 ug 0M (z + ,u;v + ,v) = 9(z + — v + ) 


(2.30) 


and defining the Wilson line following the continuous trajectory z(z + ) from y + to x + as 


U(x + ,y + , [z(z + )]) = V+ exp 


ig / dz + T • A (z + ,z(z + )) 


(2.31) 


Next, let us consider the contribution 1 to TZ k (x + , y + \ y), in which in one discrete step 
only the terms subleading with respect to ein eq. (2.29) are kept, whereas in all the 

Gift m 

other steps, only e“v is kept. This contribution reads 


Kk(x + ,y + -,y) 
N -1 


1 

iV-1 


= lim 

TV—>-oo 


' N \ ( N -1 ^ 

n d2u - n Go,k+ ( Z n+1’ u n+l) z n > u n) r 

,n=l J \ n=0 J 


n 

p= 1 l. n=p-\-l 


a n 

e n 


r-l 


/o 


f 1 , f 1 s , s bp(u p ) u o^p bp(u p ) fl - s -^2£ , ^ 

+ / ds / due n ^ s ^ + O 

do Jo N A 


fop( u p) 

JV 


p—i 

n 

n=0 


a n 

e n 


n =0 


JV -1 ,. r AT —1 

= J im / rf2u P £o,fc+(^p ,u p ;y + ,u 0 ) 7M TTe^ 

iV—>oo z ' / I 

P=1 

/*1 /*1 —s 

+ ds 
Jo JO 


"1 


i , s ^p bp(u p ) fp 

/ ds e n e n 

I o A 


, &p( u p) ^p( u p) _( s + u )2£ . ^ 

du e n e n e + ; jv + O 

N N 


bp( u p) 

N 


(2.32) 


using the identity given in eq. (2.30) as well as the relation 

J d 2 z g 0 'k+(x;z ) g 0 ,k+U',y) = 0(x + -z + ) 9(z + —y + ) Go tk +(x;y) . (2.33) 

In the continuum limit N —> +oo, the sum over p becomes an integration over z + , with 
the correspondence 


N-l 

h E - 


JV 


p=i 


an" — 


^ [ X dz+ 

y + Jy+ 


(2.34) 


Hence, the term linear in b p (u p )/N in eq. (2.32) gives a finite contribution in the continuum 

a P 

limit, whereas all the other ones converge to zero. Moreover, one has e s iv —y 1, and thus 
the integration over s becomes trivial. All in all, one gets 

i n ~i f r N - 1 

Kk(x + ,y + -,y) = Jim^ — ^ / d 2 u p £ 0jfc+ (z+, u p ; y + , u 0 ) V + < b p (u p ) \\e^ \ 

1 ^°° V p= 1 ^ l n=0 J 

= J + dz t J d 2u 1 ^o,fc+( 2: i"> u i;y + ,uo)T’ + |jd(j: + ,?/ + , [z(z+)]) igT • &Aj"(ui) j ,(2.35) 
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with the notation 


5A n (u n ) = A (z+; z n + u n ) - A (z+; z„). 


(2.36) 


Here, the V+ means that the decoration igT ■ <L4j~(ui) has to be inserted into the Wilson 
line U ( x + ,y + ; [z(z + )]) at the appropriate place in order to fulfill the x + ordering. That 
notation will be used frequently in the rest of the paper. 

These results generalise easily to the leftover contributions to the medium modification 
factor 7Zk(x + , y + ; y), with only the terms subleading with respect to e~v in eq. (2.29) kept 

dft . 

in more than one discrete step, and only e~v in the others. Adding all the contributions 
together, one obtains in the continuum limit 

^ +00 I > x + - x + 

nf(x+,y + - y) = U ab (s+, y + ; [z(s+)]) + ^ / dz f / d 4 ’ ’ ’ 

1=1 d v + dz t 

U (x + ,y + - 1 [z(z+)]) 

3 =1 

2.2.3 Expanding around the classical trajectory 

With the expression (2.37), we are now in the position to perform the semi-classical expan¬ 
sion of TZ^ b (x + , y + ; y) around the straight trajectory z(z + ) with a fixed angle associated 
with k /k + . The method is the following. First, one rewrites the finite difference 8Aj(uj), 
defined in eq. (2.36), as a Taylor series in Uj every time it appears in eq. (2.37). Second, 
one can perform the integration over each of the transverse displacements up Upon inte¬ 
gration, each power of Uj converts into a power of 1 /y/JA. This method can be pushed to 
arbitrary orders in principle, but becomes more and more cumbersome at higher orders. 
In the present paper, we limit ourselves to the order l/k +2 , corresponding to next-to-next- 
to-eikonal accuracy. The details of the calculations are given in appendix A. 

For simplicity, let us introduce the following notations: 


B\i) = igT ■ d yi A~(z + , z {z + )), (2.38) 

= igT ■ d r d y] A~(z + ,z(z + )), (2.39) 

B l3 \i) = igT ■ d yi d yj d y iA~(z + ,z(z + )), (2.40) 

B ljlm ( 3 ) = igT ■ d y id y jd y id y mA~(z + ,z(z + )) , ( 2 -41) 


IF 

3 = 1 


o,fc+l 


"A-, i-W 


V , 



where 3 = (z + , z(z + )). 

Collecting the results from Eqs. (A. 8 ), (A.ll), (A.14) and (A.17) in appendix A, one 
finds that the semiclassical expansion of the medium modification factor 7Z l jf{x + , y + ] y) 
reads 


Kk(x + ,y + ;y) = U(x + ,y + ; [z{z + )j) + i ^ X - U {1 } (x + ,y + ] [z(^ + )]) 


^ \k+)* u {4 x+ ’y + '> [*(* + )D + ° (( 


(x + -y + ) n 2 

k+ 


(2.42) 
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with the next-to-eikonal term 


^{i}(x + ,y + ;[z(z + )]) = V + U (x + ,y + -[z(z + )]) j f dzf P + _ U + \ <$* J # J ( 3 i 

I Jy+ \ x y ) 


+2 


/ X dz+ f dz+ ^_ y ^ , (2.43) 

Jy+ J z{ \ x V ) I 


and the next-to-next-to-eikonal term 


U {2} (x + ,y+;[i(z + )]) = V + U(x + ,y + -[i(z + )]) | \ f/4 ( ^+_^ ) )2 &8 lm 


+2 


pdzt [ X dz+ \ (4 p)(4 y+ ) (s^ 1 B^\ d2 )B\ Sl ) + ^6 lm B^( d2 )B lm ( dl ) 
Jy+ J zf L \ x y ) \ 


+ 


(zj-y + ) 2 

(; x + -y + ) 2 

n+ 


rx 1 px 1 

+2 / cfcj*" / / dzg" 

«/ 2^"" 2^"" 


^( 32 )^^( 3 i) + ^( 3 2)^'( 3l )) 

(4-y + )(4-y + ) 


'y 


+ 


(x + -y + ) 2 


S l W lJ (te) B\i 2 ) #( 31 ) 


+ P z i_ y l\ 2 B l (i 3) #(a 2 ) ^( 3l ) + (z 2 + y + K^| y + ) (gi (3 3) ^6^( 32 ) #( 3 i) 

{ x y ) \ x y ) ^ 

+B' l (z 3 ) #( 52 ) <y J ' z #' , ( 3 i) + 2 #^'( 33 ) #( 3a ) #'( 81 ) + 2 #( 33 ) #'( 52 ) #'( 31 )) 

+4 fdzt Pdzt p *4 r d 4\ {z *^+-Sv y+) Bi{u) Bl{53) &{n) BJ{hl) 

J ?y+ J zt J Zn j Zo L y ) 

' ' ■ . (2.44) 


+ ( 4 ^+)(z 1 +^ + ) gi(34) ^( 33 ) ^( 3l ) + ^( 3 2 ) #( 31 ) 


(x+-y +) 2 


Thanks to longitudinal boost invariance, the expansion (2.42) can be understood either as 
a large £; + expansion or as a small x + —y + expansion. The small parameter (x + —y + )/k + is 
dimensionful, but this is compensated by the appearance of two more transverse derivatives 
of the background field at each order in the expansion. 


2.3 Small angle expansion for the background propagator 

So far, we have performed the large k + expansion of TZ 1 ^(x + , y + ; y) at fixed k /k + , with 
the result given in eq. (2.42). This is enough if one is interested in hard processes i.e. 
in the Bjorken limit. However, we are mostly interested in processes in the high-energy 
(or Regge) limit. Then, we have to re-expand the expression (2.42) in the small deflection 
angle limit k/£; + —> 0. In eq. (2.42), the quantity k /k + appears through the trajectory 
z (z + ) of the (decorated) Wilson lines, defined in eq. (2.17). 

The term of zeroth order in eq. (2.42) is just the Wilson line U(x + ,y + ; [z(^ + )]) along 
the trajectory z(z + ). Its Taylor expansion with respect to \s./k + at second order is straight- 


- 9 - 


















forward to perform, and gives 


U (x + ,y + ][i{z + )}) =U{x + ,y + ]y) + (x + -y + )— U^ 01] (x + , y + \ y) 


+ ^ + -,+ ) 2 ^ U^ 2] (x + ,y + ] y) + 0 


(x+-y+)|k| N 3N 
- k+ - 9± 


where 


and 3 


% n (x + , y + ; y) = T+ U (x + , y + ; y ) f dzf ^ V j B‘ (zf, y), 
1,1 Jy+ ( x+-y +) 


(2.45) 


(2.46) 


^021 ( x+ iy + ',y) = ' p +u( x+ ,y + 'iy){ / d 4 


+ / 4- 4_ \ 2 

zi -y 


1 [ 0 , 2 ] 
+2 


x^—y^ 


& j (4,y) 


J + dz t f z+ d 4 2 (j / +ly+) 2 y ~ B *(4> y) & (4 . y)| ■ ( 2 - 47 ) 


In contrast, we need to Taylor-expand Zi{ 1 }(x + , y + ; [z(z + )]) only to first order in 
k/fc + since it appears in the first order term in eq. (2.42). According to eq. (2.43), 
£f{i}(x + ,[z(^ + )]) also depends on k /k + only through the trajectory z(z + ), which ap¬ 
pears both in the Wilson line itself and in the decorations. Hence, one gets 

U{\}{x + 1 y + ] [z(z + )]) = V + ^U{x + ,y + \y) + {x + -y + )j^ 4 [01] (x + ,y + -,y)^ 

X { f X + dz+ i ( x +- y y + +) [ Sij B 4 Z+ ^) + ( z+ ~y + )^t jl B ijl (z + ,y) 




B\ z i,y) ® l X4, y) 


+( 4 -y + ) ^ ^(4, y) B\z+, y) + (z+-y + ) ^ B\z+,y) B^(z+,y) 


+o 


(x + -y + )\k\ 

k + 


d± 


(2.48) 


Eq. (2.48) can be organised and rewritten as 


U{i}(x + ,y + \ [z(2 + )]) =^[i i0 ](^ + ,y + ;y) + ^ - k 4k4 x+,y+ ' ,y ) 


+o 


(x+ —y+)|k| N 2N 
- k+ - 


(2.49) 


3 It is important, at this point, to make a clarification on the notation of the decorated Wilson lines. 
Here after, each decorated Wilson line appears with a subscript [a, 0] where both a and /3 can take values 
from 0 to 2. In this subscript a stands for the order of expansion around the classical trajectory and p 
stands for the order of small angle expansion. For example, W( 0 ,, is the decorated Wilson line originating 
from the zeroth order expansion around the classical trajectory and the first order small angle expansion. 
We also would like to emphasise that Uk -n and W[i,o] were called and U( 2 ) respectively in Ref. [1], 
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where 


U[ifl\{x + ,y + ]y) = 'P+ U(x + ,y + ; y)| J^ + dz f ^(^y) 

+2 [ X d zt [ Zl dzt {Z i~ y+ + \ B\zt,y) B\z+, y) 

■V V (x+-y+) 


(2.50) 


and 


w [\,i]( a:+ )y + ;y) = ^(^ + ,y + ;y)s / ^1" 


x+ — 


& il (zt,y) 


+2 / dzj* - / ciz^ 


-y 


+ 


x^-y 

(zj-y + )( z 2-y + ) 

(x + -y + ) 2 


B j (z+, y) B^+y) 
^(4,y) B\z\ ,y) 


+^fi i (4»y) ^^(4,y) + ^^ 7 (4,y) B^z+y) 


+2 / dzi 1 " / dz^ / dz^“ 


(23+ (J / + + -y+) 2j/+) gt( " 3+,y) 6 ' (4,y) ^ ( ^’ y) 


+ (Z2+ (, x + + i ( j()2 ?/+) gJ (4 1 y) (^(4.y)^y) + ^(4,y)s ! (h + j) 


. (2.51) 


For the last term in eq. (2.42), which is already of next-to-next-to-eikonal order, we can 
drop completely the k /k + dependence, or equivalently make the replacement z(z + ) t —> y. 
We thus write 


where 


U{ 2 }(x + ,y + -,[z(z + )]) =U[ 2 ,o](x + ,y + -,y) + O ^, 


W[2,o] (^ + , y + ; y) = W {2 } (.x + , y + ; [y]), 


(2.52) 


(2.53) 


which can be read from eq. (2.44). 

All in all, after both the semi-classical and the small angle expansions, the medium 
modification factor TZ^ b (x + , y + ; y) for the propagator reads, at next-to-next-to-eikonal ac¬ 
curacy, 

Kk{x + ,y + -,y) =U(x + ,y + ] y) + ^ — U{ 0A] (x + ,y + ; y) 

.(x+-y + ) + + ^ , (x+-y + ) 2 k*k? {j + + 

+ * 2k+ U [ i, 0] (x + ,y + ;y) +-^- U^ 2] (x + ,y + -, y) 

{x + —y + ) 2 k* , , , (x + — y + ) 2 I , 

+* 2(fc+) 2 w [i,i](* + .y + ;y) - 4(fc + )2 %o](^ + ,y + ;y) 


+o 


(x + -y+) |k| 
k + 


di 


+ 0 




(2.54) 
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where the various decorated Wilson lines are defined in Eqs. (2.46), (2.47), (2.50), (2.51), 
and (2.53). 

3 Single inclusive gluon production in pA collisions 

The eikonal expansion performed in the previous section at the level of the gluon back¬ 
ground propagator can be applied to single inclusive gluon production in pA collisions at 
high energy. This process can be treated in the CGC effective theory assuming that the 
gluon is produced at mid-rapidity. 

Within the CGC framework, following the same idea developed in section 2, we describe 
the highly boosted left-moving nucleus by the background field given in eq. (2.4) with a 
finite longitudinal support from x + = 0 to x + = L + . The right-moving proton, on the 
other hand, is considered to be dilute and described by a classical color current 

f{x) = ^~jt{x) (3.1) 


which is localized at x~ = 0. 

Let us consider a proton-nucleus collision with some impact parameter B, taking the 
centre of the nucleus as the reference point for the transverse plane. Then an arbitrary 
point x in the transverse plane is at a distance |x — Bj from the center of the proton and 
at a distance |x| from the center of the nucleus. In this setup, the classical color current 
j + (x) is written as 

ji{x) = 8{x~)U ab {x + ,-oo-yi) p b (x —B), (3.2) 

where p b is the transverse color charge density inside the proton before the interaction with 
the nucleus, and U ab (x + , — oo; x) is the Wilson line implementing the color precession of 
these color charges in the background field A~(x + ,x) of the nucleus. 

The single inclusive gluon production cross-section for a pA collision is given by 

E (<IUS(fcB)|*) ) , (3.3) 

A phys. ' ' A 

where A, a and k = (fc + ,k) are the polarization, color and momentum of the produced 
gluon. At leading order in the coupling g , the gluon production amplitude, is given 

by the LSZ-type reduction formula 

M a x (k, B) = ef J d A x e ik x O x A“(x) , (3.4) 

where A“(x) is the retarded classical field [40]. Since we are interested in dilute-dense 
scattering, the background field A~(x + ,x.) that is describing the dense nucleus is 0(l/g) 
and the color charge current (x) that is describing the dilute proton is 0(g). The retarded 
classical field A“(x) appearing in the reduction formula contains both the background field 
A”(x + ,x) and small perturbations on top of it, generated by the color current j+(x). 
However, since in the light cone gauge A + = 0 the physical polarization vectors satisfy 
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e A* = 0) only the transverse components of the field perturbation due to j+( x) contributes 
to eq. (3.4), which can be written at leading order in g as 

M a x {k,B,) = e l x lim e lk ~ x+ [ d 2 x e"* k x [ d A y Gl+{x;y) ab j£(y). (3.5) 

a;+->+oo J J 

The gluon production amplitude .M“(fc, B) can equivalently be written in momentum 
space as 

M a x (k, B) = f ^ e _iq ' B M?(k, q) p b ( q) (3.6) 

by defining a gluon-nucleus reduced amplitude M°\ {k, q), and p( q) through the Fourier 
transform 

p“(y - B ) = / P°(q) • P- 4 * * 7 ) 

It is clear from eq. (3.3) and eq. (3.6) that one needs to average over the color charge 
densities at equal momenta, (p a (q)p b (q)) p , in the calculation of the single inclusive gluon 
production cross-section, due to the integration over the impact parameter B. It has been 
shown and discussed in detail in Ref. [1] that this correlator is indeed related to the 
unintegrated gluon distribution (^(q ) 1 in the projectile via 

(p a ( q)* P b ( q)) p = -J?—[ q2 Vp{q) • (3.8) 

3.1 Corrections to the eikonal limit at the amplitude level 

The gluon-nucleus reduced amplitude M. ab {k, q) can be decomposed into three contribu¬ 
tions as 


Mf(k, q) = M a b b ef:X (k, q) + Mi b jX (k, q) + q), (3.9) 

in which the gluon is radiated by the color current respectively before, during or after 
the interaction with the background field A ~ (x + , x) . Each contribution can be written 

4 Here, we indicate only the dependence of <p p on the transverse momentum of the gluon, both for 

simplicity and by consistency with the fact that our calculation is only a leading order one in the coupling. 

On the one hand, due to the boost-invariance of the classical Weizsacker-Williams gluon field sourced by 
the current ip p ( q) does not depend on the k + of the gluon at this order. On the other hand, when 

including radiative corrections, y p (q) would acquire dependence on a factorization scale regulating the soft 
divergence. In order to resum optimally the small-a; leading logs for the projectile, one should then take 
for the factorization scale a value related to the k + of the gluon. This situation is analog to the case of 
collinear factorization, where the parton densities acquire a dependence on the hard scale through the choice 
of renormalization scale, which is absent in a leading order calculation in the parton model. 
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explicitly as 


M? ef , x (k,ci) = efe ik - L+ i J d 2 y e™ (-2) ^ f d 2 ze~ ikz g$ (L+, z; 0, y), (3.10) 


= £\e lk L+ i J d 2 ye iq ' y ^| ciy + |a yl J d 2 ze lk z 0^(L + , z;y) 

xM c V-0;y), (3.ii) 


q) = 4* e ik ' L+ i j d 2 y 2 ^ ^ a6 (L + , 0, y) . (3.12) 

The ” aft ” contribution to the total gluon-nucleus reduced amplitude does not involve 
the background propagator Gi+(x,y), thus it can be kept as it is written in eq. (3.12). On 
the other hand, the ”6e/” and ”m” contributions to the total reduced amplitude involves 
the background propagator whose eikonal expansion has been performed in section 2. Thus, 
one can use the expanded expression of the background propagator G%+(x , y) in eqs.(3.10) 
and (3.11) in order to calculate the corrections to the eikonal limit of the ”6e/” and ”in” 
contributions to the total gluon-nucleus reduced amplitude. 

Let us first consider the "bef" contribution. By using eq. (2.16), M be j X (k. q) can be 
written in terms of the medium modification factor lZk(x + : y + ;y) can be written as 


q) = el* J d 2 y (-2) ^ Uf (L+, 0; y). (3.13) 

One can now use the expression of the expanded IZ'jf (L + ,0',y) that is given in eq. (2.54) 
in eq. (3.13) to write the ”bef ” contribution to the total gluon-nucleus reduced amplitude 
at next-to-next-to-eikonal accuracy as 

<Afeq) = ijd 2 y (-2) |ZL(L+, 0;y) + ^ k* U( 0>1] (L+, 0; y) 

+* ^[i,o](^ + ,0;y) + kV W” 2] (L+ 0;y) + % - (^j k^ U 3 {11] {L+, 0; y) 

-j(j+) M[ 2 , 0 l(i + .° ; y)} . (3.14) 

The remaining contribution that involves the background propagator G%+ (x, y ) is the 
”m” term M^ x (k, q), given in eq. (3.11). Inserting eq. (2.16), into eq. (3.11), one gets 


Tab 


l 


r L+ 


q) = 4* * d 2 y e l(k q) ' y ^T / 




(<9 y ,: -^ kl )^r( L+ >2/ + ;y) 


xwV.°;y)' 


(3.15) 


The presence of an explicit factor of l/fc + in eq. (3.15) suggests that it is enough to 
keep the terms of O (1) and of O (L + /k + ) in the expansion of the medium modification 
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factor TZ < f c b (L + 1 0;y) in order to get the expression of the M^ nX (k, q) at next-to-next-to- 
eikonal accuracy. Keeping this in mind and using the expanded expression of the medium 
modification factor 7£^ b (L + ,0;y) given in eq. (2.54), M inX (k,q) can be written as 


q) = 4* i / d 2 y ^ dy+ e^ +fe “ (d y , - iW) 


1 




+(L + -y + )^W [0A] (L + ,y+-,y) + ^l T (L + -y + )U [lfi] (L + ,y + ;y) 


k + jo 
2k + 1 


U(L + ,y + ; y) 


ZY c6 (y + ,0;y). (3.16) 


Note that eq. (3.16) contains terms with the transverse derivative acting on both the 
usual Wilson line U ab (L + , y + \y) and also on the decorated Wilson lines U^,(L + , y + ; y) 
and (L + ,y + ] y). In order to have an explicit and complete expression for the ” in ” 
contribution to the gluon-nucleus reduced amplitude at next-to-next-to-eikonal accuracy, 
one should calculate these terms explicitly. Transverse derivative of the usual Wilson 
U ab (L + , y + ; y) line has been calculated in Ref. [1] and the result reads 

d r U ab (L + ,y + - y) = f dz + [U{L + ,z + ] y) (igT ■ d r A~(z + , y)) U(z + , 0; y )] ab .(3.17) 
Jy+ 

On the other hand, the transverse derivative of the decorated Wilson lines have not ap¬ 
peared at next-to-eikonal accuracy, thus they have not been calculated in [1] . Given the fact 
that the decorated Wilson lines are composed of the usual Wilson lines and field insertions 
between these usual Wilson lines, it is straightforward to calculate the transverse derivative 
of the decorated Wilson lines. When the derivative is acting on the field insertion, it simply 
increases the number of derivatives acting on the background field Ar(z + ,y). When it is 
acting on the usual Wilson line, it introduces an extra field insertion as shown in eq. (3.17). 
Then, one obtains 

d y iU [o,i ]( L+ >2/ + ;y) = 'P+ld(L + ,y + ; y) dz + _^ + B l] (z + ,y) 

+ V+U{L + ,y + ] y) f dw + B l {w + , y)f dz + ~ V B 3 (z + , y). (3.18) 
J y + J y + L+-y+ 

Here we only present the transverse derivative of one of the decorated Wilson lines but 
the same arguments hold for the other one and it can also be calculated straightforwardly. 
All in all, the ” in ” contribution to the gluon-nucleus reduced amplitude M in x (k, q) at 
next-to-next-to-eikonal accuracy can be written as 


M? nyX (k, q ) = e\*i I d 2 y e * {q - k) y| 2 p(l- e ^ +fe -)z7“ b (L+,0;y) 


L+ 


L+V \i 


+^T W[ 0 ?](^ + .°;y) + ( ) 4 (k 2 ^- 2 k l k^)^ A) (L+, 0 ;y) 


0; y) + 4^( C )( l+ ’ 0; y) 


T ab 


(3.19) 


- 15 - 









Here, Zd^(L + , 0;y), U^(L + , 0;y) and U l ^(L + , 0;y) are the new decorated Wilson line 
structures that are defined as follows: 


r L+ 


f Z + \ 


Z4$(L+ 0;y) =V + U ab (L + ,0-y) J dz + B i (z + ,y), 


(3.20) 


,L+ 


+ \ 2 


U^ b {L+, 0; y) = V+ U ab {L + , 0; y) <j \d ij 5 lm + 8 il 6 jm + S im 8 jl J j dz + ['—) B lm (z + ,y) 


+2 


rL+ rL+ 


6 ij d lm + § il 5 jm + 5 im 6 jl^ j ^ + <fe+ ) B^+y) S m (4,y) 

fL+ /-L+ 


+4 J Q dzt f d4 Zl {Z { l + y V ] &{4,y) B\zf, y)}, 


(3.21) 


^(C) b (-^ + : 0; y) = P+ U ab (L + , 0; y) / dz + & l & l {z+, y) 


/■L+ /-L+ 

+ / dw + / cfe + 
Jo Jo 


-(- \ 2 _|_ 

6>(u + -z + ) ( ) +^ + -u + ) 7 ^(2z + -r + ) 


L+ J ' ~ K ~ ~ ' (L+) 2 

x y) B>+ y) 

+2 ^dzt fjdzi \B^(z+,y)B\zt,y)+B^z+,y)B^(zt,y) 


L+ J 

rL+ f L+ rz+ 

+2 / du + / dz.t' / dzi* - 

Jo Jo Jo 


% + -4)(ff) +^i + -^ + )^2 (2^i + -^ + ) 


xB l (^y)B 3 (4.y)^(4.y) • 


(3.22) 


Finally, by using eqs. (3.12), (3.14) and (3.19) the total gluon-nucleus reduced ampli¬ 


tude A4x (k] q) at next-next-to-eikonal accuracy can be written as 

Mf (t; q) = ie\* J d 2 y { 2 (p - U(L+, 0; y) 


+ 


+ 


L+ 

k+ 

L+ 

k + 


q*k J 


jy - 2 ^-)"[o,.]( i+ ' 0;y) - ^ "[i.o](r + .0;y) 

- a^k-V ttf 0 ', 2| (i + ,0;y) -i^-k :C' l l| (7.\li;yi + ( '* M[2,0](i + ,0;y) 


q 


q 


2 q 2 




+ 4 (k 2 <W - 2k*k- 7 ) W^(L+,0;y) + — W^(L+,0;y) + - ZJ* C) (L+, 0; y) 


a6 


. (3.23) 


3.2 Gluon production cross section beyond eikonal accuracy 

The single inclusive gluon production cross section and light-front helicity asymmetry of 
the produced gluon are the two observables that have been discussed in detail in Ref. [1] 
at next-to-eikonal accuracy. For the single inclusive gluon production cross section, it has 
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been shown that the next-to-eikonal terms vanish and the well known result - the strict 
eikonal limit of the A:j_-factorized formula - has been obtained. On the other hand, the strict 
eikonal terms for the light-front helicity asymmetry have been shown to vanish, leaving the 
next-to-eikonal contributions as the leading terms for this particular observable. 

We have calculated the gluon-nucleus reduced amplitude M x (k: q) at next-to-next- 
to-eikonal accuracy in the previous section. In this section, our aim is to use this result to 
calculate these two observables at next-to-next-to-eikonal accuracy. 

The single inclusive gluon cross-section, eq. (3.3), can be written in terms of the 
gluon-nucleus reduced amplitude M x (k; q) as 

(2,r>3 {2k+) = / (W l)’^ (q >), E ("Ilf q )) A ■ (3.24) 


By using the relation, given in eq. (3.8), between the correlator of two charge densities 
in the projectile (/O c (q)*/5 6 (q)) p and the unintegrated gluon distribution of the projectile 
<^p(q), the single inclusive gluon cross-section can be written as 


,+ da 
dk+ d 2 k 



1 

N?-l 


X]( 7W A b feq) t - M A 6 (^q)) A • 


(3.25) 


The second observable that we consider is the light-front helicity asymmetry of the 
produced gluon. The calculation of this asymmetry is almost identical to single inclusive 
gluon production cross section, except that instead of summing over the helicity A ± 1 of 
the produced gluon, one takes the difference between the A = +1 and A = — 1 contributions 
i.e. 


, da + _ . da 

dk+d 2 k dk+d 2 k 



X ~4~ jv2_i A (• M A > fe c O t •M?& c 0) a - (3-26) 

c A phys. 


In the calculation of both single inclusive gluon cross-section eq. (3.25) and light-front 
helicity asymmetry eq. (3.26), one needs to take the square of the gluon-nucleus reduced 
amplitude Ai x (k, q). At the squared amplitude level, we keep the A dependence explicit 
so that we can apply the result to both observables in eqs. (3.25) and (3.26). Keeping this 
in mind, the square of the amplitude can be expanded as 


M x (/c,q) f M.f (k, q) 


Sh 

+ M a x \k, q YM a x (k, q) 


+ M° x {k, q ) ] M x (k, q) 

+ • • • , 

NNE 


NE 


(3.27) 


where the terms on the right hand side of eq. (3.27) stand for the eikonal terms, next-to- 
eikonal terms and next-to-next-to-eikonal terms respectively, and the dots for higher order 
terms. 

By using the expression for the gluon-nucleus reduced amplitude M. x (k. q) given in eq. 
(3.23) each contribution to the squared amplitude can be written explicitly. The eikonal 
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contribution reads 


M?(k,^M a x b (k, q ) ^ = 4*4 / d 2 y I d 2 y ' e i(y-y')-(q-k) 


jab 


x4 C*(k, q) C J (k,q)tr [^(yO^y) 


(3.28) 


Note that we have dropped the longitudinal coordinate dependence of the Wilson lines for 
convenience but one should keep in mind that these Wilson lines run from 0 to L + in the 
x + direction i.e. U{ y) = U(L + , 0;y). We will use this notation for the rest of the paper. 
Moreover, we have introduced the shorthand notation C*(k, q) for the coefficient of the 
Wilson lines which is defined as 


c‘(k. q) -(*-*). 

Similarly, the next-to-eikonal contribution reads 

' ab '' ab ' = 4*4 [d 2 y / , d 2 y , e <(y_y ' ),(q “ k) 2 


(3.29) 


Ml (k, q) f Af“ (k, q) 

x C J (k,q) C ll (k, q) tr \lA ] (y')U l [0 ^(y) 


L+ 

k+ 


+ C l (k, q) C h (k, q) tr U (y )U 1 ^ (y 7 ) 


-ic 3 { k,q) ^ tr ^ t (y , ) w [i,o](y) +*c*( k ,q) ^ tr u (y) u [i,o](y') 


(3.30) 


where C l ( k, q) is defined in eq. (3.29) and C t3 { k, q) is defined as 


C ij (k, q) = ( ^ -2k ?: ^ ) . 


(3.31) 
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Finally, the next-to-next-to-eikonal contribution can be written as 

M?(A:,q)t^?(fe,q)| NNE = 4*4 J £ y j £ y > ^(y-y'Mq-k) 2 

x |c J (k, q) - 2 ^ k k m tr [^(yOwjo^y) “ * ^2 k * tr [^(yOwfi^Cy) 

+ 1 ^2 tr wt (y / )W[ 2 , 0 ](y) + 1 (k 2 <5* J - 2 k*k 7 ) tr [w t (y / )^( A) (y) 

+ j tr [w t (y')^(B)(y) + I tr [w t (y / )W ( * c) (y) 

+ C*(k, q) - 2 ^ k ^k m tr W(y)W[j* 2 J (y') +i^k ; tr ^(y^^y') 

+J ^ 2 tr w (y) w [2,o](y / ) _ ^ ( k2( ^ - 2kJk? ) tr w (y) w (l)(y / ) 

+ ^ tr U(y)U^(y') - J tr U^U^y’) (3.32) 

+ ^ C h (k,q)C mj (k,q) tr ^ ^(y^'^y') + J ^ tr W[ 1 , 0] (y)W [ t 10] (y / ) 
-^C^( k ,q)tr W[i ) o](y)Wj 0 | 1] (y / ) tr Wf 0)1] (y)W [ t li0] (y / ) |. 


The next step is to perform the averaging (■ ■ - )a over the target fields. Let us introduce 
the variables r = y — y* and b = ^(y + y ; ) and define the usual adjoint dipole 


0( r) = / d 2 b 


= / d 2 b 


iV c 2 -l 


N? - 1 


*K b +iy( b -i)]> 


the adjoint dipoles with one decorated Wilson line Ul "i 


2J\/a’ 


(3.33) 


M (r > = f d2b vyrr ( tr Ha ( b+ D ( b - D] )a 

= j d2b jvjvi ( tr [<« ( b + 1 ) u ( b -1)]) 


(3.34) 


and analog operators for the case of the decorated Wilson lines Z^ A )> M(b) an d ^(C) - Finally, 
let us also define the decorated adjoint dipoles formed from two decorated Wilson lines, as 

=J d " b Nf ri (*■■ Hk ( b +D «w?' (»- 1 )} ), 

= f^ b W^{ tI HA( b+ l) , 4 iA{ b -l)]) A - < 3 ' 35 > 
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Here, ■ - )a stands for the averaging over the target fields. From the cyclicity of the trace 
in eq. (3.35), one obtains the identity 




(r) 




(- r )- 


(3.36) 


Now, by using the definitions of the usual and the decorated dipole operators eqs. 
(3.33), (3.34) and (3.35), we can write the square of the gluon-nucleus reduced amplitude 
already averaged over the target fields separately as eikonal and beyond eikonal contribu¬ 
tions. The eikonal contribution simply reads 


/ -W?feq) t M?feq) 


m - 1 


= 4*4 [ d 2 r e ir ' (q ” k) 4C*(k, q) C j ( k, q) 0( r). (3.37) 

F, J 


Similarly, the next-to-eikonal contribution can be written as 


(fc,q) f 7Wf (fc,q) 


N} - 1 


= 4*4 f d 2 re-( q - k )2^^C4k,q) 


NE 


k + 


(3.38) 


q 


C‘*(k,q) q Oil] (r)-i^0 [liO] (r) 


q 


+ C*(k, q) 


C lj ( k,q) O l m (-r)+i^O m (-r) 


q 


Finally, the next-to-next-to-eikonal contribution reads 


q) 


N 2 - 1 


= 4*4 / rf 2 re ir (q_k) 2 (SY 

NNE 4 V K J 


X C J (k,q) 


- 2 i k‘k”*0|” ,(r) - i | k 1 Oj u] (r) + i O |2 , 0] (r) 


+ 2 (kV-2k-k‘) Oj A) (r) + -0;' m (r) + 2o; cl (r) 


+ C‘(k,q) 


- 2 51 k'k” Oj™ 2| (-r) + i 3) k' 0[ L1| (-r) + O|2 , 0| (-r) 

(kV - 2k J k ; ) Oj A) (-r) + | o* j(-r) - J 0 2 c) (-r) 


1 

+ 2 


c“(k,q) c™ 2 (k.q) O|-"; |i|0il| (r) + 0,i,„ ]lM (r) 


q' q y 


* q 2 ^ q) ^[ 0 , 1 ];[ 1 , 0 ]( r ) * q2 ^ q) C) [o,l];[l, 0 ]( r ) 


• q' Mi, 


(3.39) 


It is possible to further simplify eqs. (3.38) and (3.39) by considering symmetry prop¬ 
erties of the decorated dipole operators. The presence of the gluon background field A~ 
and of the projectile current )+ explicitly breaks Lorentz invariance. But for the decorated 
dipole operators, which are independent of the projectile current, the averaging 
over the background held partially restores Lorentz symmetry. In particular, it restores 
the symmetry under rotations in the transverse plane, around the center of the target. 

Hence, the decorated dipole operators should be covariant under such rotations. More 
precisely, operators like 0[i jO ](r) or O[i,o] ; [i,o]( r ): which have no transverse indices, should 
behave as scalars, whereas operators with one transverse index like O^ 0 ^ (r) should behave 
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as vectors, and the ones with more indices should behave as higher rank tensors. More¬ 
over, since the b integration has already been taken in eqs. (3.33), (3.34) and (3.35), the 
decorated dipole operators depend only on one transverse vector, r. Thus, under the trans¬ 
formation r —> — r (which is a particular rotation in the transverse plane), the decorated 
dipole operators with an even number of transverse indices are invariant, whereas the ones 
with an odd number of indices flip sign. 

By using these symmetry properties the next-to-eikonal, eq. (3.38), and next-to-next- 
to-eikonal, eq. (3.39), contributions can be arranged as 




N? - 1 


= 4*4 f d 2 r e * r '(q- k ) 2 ^ 


NE 


C4k,q)C4k,q)-C’(k,q)C4k,q) 




k+ 


C4k,q)4-C l (k,q)^ 

q- q- 


(3.40) 


O[i l0] (i 


and 


f M a x {k,q^M x (k,q) 


N? - 1 


NNE 


= 4*4 / d?v ( Tf) 


c j ik<,j y-c'fk.<!)',, 

q- q- 


- 4 k l k m Oj’™ (r) - 2i k l Oj ltl] (r) + O m (r) 


+ 2 

1 

+ 2 


C j { k, q) (k 2 5 il - 2k l k l ) + C( k, q) ( k 2 d jl - 2k>k l 


0{ A) (r) 


C j ( k,q) S il +C i (k,q)6 jl 


k m O%(r) + iO\ c) (r) 


+ 


+i 


C L \ k,q) C-4k,q) Ofr 1]M (r) + ^ ^ 0 Mm (r) 


'(B) 


q* q J 


(C)V 


^^(k,q) + ^C4k,q) O 


[ 0 . 1 ];[ 1 , 0 ] 


(r) • 


(3.41) 


The final expressions for the next-to-eikonal, eq. (3.40), and next-to-next-to-eikonal, 
eq. (3.41) contributions at the squared amplitude level are ready to be substituted to 
the aforementioned observables. Before we proceed further, we would like to point out an 
important realisation about the transverse momentum structures of the eikonal, next-to- 
eikonal and next-to-next-to eikonal terms. 

The transverse momentum structure of the strict eikonal and next-to-next-to-eikonal 
terms are symmetric under the exchange of i j unlike the next-to-eikonal terms whose 
transverse momentum structure is anti-symmetric under the same exchange. This affects 
directly the observables we are interested in. 

In the calculation of the single inclusive gluon cross section one should sum over the 
gluon polarizations, which leads to 


E4*4 = ««. (3.42) 

A 


Since the transverse momentum structure of the next-to-eikonal terms are antisymmetric 
under the exchange of i j, this contribution to the single inclusive gluon cross-section 
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vanishes leading to the following result at next-to-next-to-eikonal accuracy 


k + 


= ^ (q> T /* e "' <q ' k ’ 4C ‘ (k ’ q) C ‘ (k ' q >° (r) 


+ (w> 


2^C l (k,q) 


4 k z k m 0j’™ (r) - 2i k* 0| ltl] (r) + 0 [2 , o] (r) 


+C*(k, q) 

+C K (k,q; 


*(k <y* - 2 k*k‘) Oj A) (r) + k m 0 (b)( r) + i 0( C) ( r ) 


C"''(k, q) 0[o,i] ; [ O ,i]( r ) + 2 * 2 C ’[o,i];[i,o]( r ) 


+ ^2 °[l,0];[l,0]( r ) 


.(3.43) 


On the other hand, the light-front helicity asymmetry, as mentioned before, is calcu¬ 
lated by taking the difference between the A = +1 and A = —1 contributions. In this case, 
one uses 


£ A£"4 = i £ « (3.44) 

A 


where e*- 7 is the antisymmetric matrix with e 12 = +1. Thus, the light-front helicity asym¬ 
metry gets the contribution from next-to-eikonal terms. The strict eikonal and next-to- 
next-to-eikonal terms vanish due to the fact that the transverse momentum structure for 
these terms are symmetric under the exchange of i j. The final expression of this 
asymmetry reads 


k + 


da + 

dk+ d? k 

L+ \ 


- k + 


da 


dk+ d? k 


d 2 q 

( 2 ^ 


M q) q 2 / ^ e ' r ' (q_k) 


— j e lJ C 7 (k, q) U C Ll (k, q) 0[ o>1] (r) + ^ 0 [1)O] (r) 


(3.45) 


4 Conclusions 


The retarded gluon propagator in a classical background field is one of the most important 
building blocks of the high energy dilute-dense scattering processes and also of medium- 
induced gluon radiation. Corrections beyond the eikonal approximation for this propagator 
were first considered in Ref. [1]. The corrections associated with the finite length of the 
medium were calculated at next-to-eikonal accuracy, meaning at order O (If^ 2 )- These 
corrections involve new operators referred to as decorated Wilson lines, as they include 
insertions of gradients of the background field along the path. 

In this paper, we extend our study of finite width effects on the retarded gluon prop¬ 
agator, including next-to-next-to-eikonal corrections that are order (iFl) ' The new 
operators that appear at next-to-next-to-eikonal accuracy are also decorated Wilson lines 
but with higher number of insertions of the gradients of the background field or with higher 
derivatives. 

The eikonal expansion performed at the level of the gluon background propagator is 
then applied to high energy dilute-dense scattering processes within the CGC framework. 
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Two different observables have been analysed, in pA collisions at midrapidity, within this 
framework: the single inclusive gluon production cross section and the light-front helicity 
asymmetry of produced gluons. For the single inclusive gluon cross section, it has been 
shown that the next-to-eikonal terms vanish and the first non-vanishing corrections to the 
strict eikonal limit that appear at next-to-next-to-eikonal order have been calculated. On 
the other hand, for the light-front helicity asymmetry, it has been shown that both the 
strict eikonal terms and next-to-next-to-eikonal terms vanish and the leading contribution 
to this observable turns out to be the next-to-eikonal terms. This result shows the analogy 
between the twist expansion of the hard processes and the eikonal expansion of the high- 
energy processes. 

The decorated dipole operators appearing in both of the observables at next-to-eikonal 
and next-to-next-eikonal accuracy are expected to have rapidity divergences. Therefore, 
understanding the low-x evolution of these operators is a complementary extension of the 
analysis of the CGC beyond eikonal accuracy. This issue is left for future studies. 

The eikonal expansion that is considered in detail in this paper has further applications. 
It can be applied to other high energy dilute-dense processes such as DIS and single inclusive 
gluon production in the hybrid formalism, or to jet quenching physics. Such applications 
are also left for future studies. 
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A Details of the semi-classical expansion 

In this appendix, we calculate the semi-classical expansion of lZ^ b (x + , y + ; y) up to order 
l/k +2 starting from the expression (2.37), using the method outlined in section 2.2.3. For 
simplicity, we consider separately each term with a given index l from the series represen¬ 
tation (2.37) of 77j) fe (x + , y + ; y). Note that terms with l > 5 can only give contributions 
more suppressed than \/k +2 in the k + —> +oo limit. 
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A.l l = 1 term in eq. (2.37) 

The term with / = 1 in eq. (2.37) can be written as 


/+ dz+ j d2u ^o,fc+(2 + ,u;y + ,0) 

xP+jw (x + ,y + ; [z(z+)]) ig T • [A~{z + , z(z + ) + u) - M~(z + , z(z+))] | 

= V + U (x + ,y + ;[z(z + )j) J^ dz + ig T ■ J d 2 u g 0)k+ (z + , u; y + , 0) (A.l) 


’y~' 

u' IjJ iTxTu* u*u J u i u m 

u'a, + —d t d, + — r d l d 3 d l + --- didjdidm 


A (z + , Z(z + )) \ . 


Note that in this equation, each partial derivative is implicitly with respect to y and they 
are acting only on the background field A ~. Also note that the free scalar propagator Go,k+ 
is even in u, so that any odd power of u appearing in the expansion vanishes. Moreover, 
due to azimuthal symmetry, we can make the replacements 


, , 2 

u U J — u , 

2 

lu 2 l 2 

uWu m ^ 1— L^S lm + S il S jm + 5 im 5 jl ). 

8 

Thus, the l = 1 term in 1Z% b (x + , y + ; y) can be written as 
Kf (x + ,y + ;y)i=i = V+ U (x + ,y + ; [z(z + )]) igT • f dz + 

Jy+ 

1 [ d\A~ {z + A{z + ))\ [ d 2 u (u 2 ) Q ok +(z + ,u;y + ,0) 


(A.2) 

(A.3) 


X< 4 

1 

+ 64 


(■dlfA (z + ,z(z + )) / d z u (u z )' z g oM (z + ,u-,y + ,0) 


2\2 


(A.4) 


Using the expression (2.14) for the retarded free scalar propagator, it is straightforward to 
obtain the following relations : 


I d 2 u g o k+ (z + ,u-,y + ,0) (u 2 ) = 2 i > Z 
d 2 ug oM (z + ,u-y + ,0) (u 2 ) 2 = -8 


k+ 

(z + - y+) 2 

(. k +) 2 


(A.5) 

(A.6) 


Then, by using eqs. (A.5) and (A.6), the l = 1 contribution to 1Z% b (x + , y + \ y) reads 
/ Rf(x + ,y + -y)i =1 = V+ U (x + ,y + ; [z(z + )]) igT- j ^ dz + j V" - d 2 ± A~ (z + , Z(z + )) 


[z^ - y 

8(k+y 


+ ^2 r 


(dl) 2 A~{z + A(z + )) 


(A.7) 


- 24 - 




















or, using the notations introduced in Eqs. (2.39) and (2.41), 

K a k(x + , y + - y) i= i =V+U (x + , y + ; [z(z + )]) 

x J y l “ %k+y 2 sij6lm Bijlm(3l) 


(A.8) 


A.2 l = 2 term in eq. (2.37) 

The term with l = 2 in eq. (2.37) reads 

/ x + PX~^~ P P 

dzf j dz£ J d 2 ui j d 2 u 2 G 0ik +(z^,u 2 -,zf,u 1 ) 

X V + U (x + ,y + -,[z(z + )}) igT-5A 2 {\i 2 ) igT • 8A± (ui) G 0 ,k+(zf ,u 1 -,y + ,0) 

/X~^~ PX~^~ P P 

dz+ J dz£ J d 2 ui J d 2 rG oM (zt,r-,zt,0) (A.9) 

x V + U (x + ,y + -[z(z + )}) igT • 5A 2 (ui + r) igT-5A[( ui) C/ 0jfc+ ( 2 ^, ui; y + , 0) , 

where r = 112 — ui. Taylor expanding the (5A~s with respect to ui and r, and dropping 
the terms vanishing due to azimuthal symmetry, eq. (A.9) becomes 


Kt {x + ,y + ]y)i =2 
= V+U (x + ,y + ; [z(^ H 

x Go,k+( z i’ u Tiy + ’°)\ u i & l (i 2 ) u i ^ J (3i) + 
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1 

+ 2 


dz i I dz 2 l d2u 1 /d 2 * V(4. r ; z n°) 

^(32) Uf B m (u) 
(A.10) 


— u) Ui u( + 1 r'rhi 
6 1 1 1 2 
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uXi + rV B lJ ( 32 ) ^u'u T B lm ( 5l ) + u\ B\i 2 ) 


J,, 


6 



Then, using eqs. (A.2) and (A.3) and integrating over ui and r by eqs. (2.30), (A.5) and 
(A.6), we get 


Kk , {x + ,y + \y)i= 2 = V+ 


U (x + ,y + ][i(z + )}) j dzf j ^ &(i 2 ) fffoi) 

J y~^~ «/ Z -^ v 


(4-y + )( z t~y + ) . 1 


2 (k+y 

(z+-y+) 2 r 


s jl B { i\ 32 ) ^( 31 ) + 5 lm B lm (n) 


2 {k+f 


B i (32)^B ij7 (3i)+fi <j (3 2 )B ij '(3i) 


(A.ll) 
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A.3 l = 3 term in eq. (2.37) 

For 1 = 3, one has 

K-f(x + ,y + ;y)i =3 = V+ U (x + ,y + ; [z(z + )]) f dz\ f dz£ I dz£ 

J y+ J z+ J z~2 

X / l / d2U3 ^°> fc+ ( Z 3" ,U 3 ; ^ ,U2 ) £°,fc+(^ U 2 ;Z l’’ U l) 

x g oM (z+,ur,y + ,0) igT ■ 5A% (u 3 ) igT ■ 5A^ (u 2 ) igT-6A^{ Ui) 

= T + (x + ,y + ; [z(z+)]) [ dzf [ dz£ [ dz£ 


z + 

Z 1 


x J d 2 s J d 2 r J d 2 ug 0tk +(z£,s-,z£,0) G o k +(z^,r; zf ,0) G oM {z^,u-y + , 0) 
x igT • A 4.3 (s + r + u) igT • 8A2 (r + u) igT ■ 5A"[ (u), (A. 12) 

where s = U 3 — U 2 , r = U 2 — ui and u = ui. Taylor expanding with respect to s, r and u, 
and dropping the terms integrating to zero, one gets 

K% b (x + ,y + ;y)l=3 = p+ U (x + ,y + ] [z(z+)]) f dzf [ dz£ [ dz£ 

Jy+ Jz+ Jz+ 

X J d 2 s J d 2 r J d 2 uG oM (z£,s;z£,0) G 0:k +(z£,r;zl,0) G 0)k +(zf,u-,y + ,0) 


x < -u'nVu"' 


^(33) B \ 32) B m (i 1) + ^(33) B*\ 32) B' m { 31) + ^(33) #(32) ^ m ( 3 i) 


IrVu'u” 


B il (, 33) ^(32) + ^(33) ^(32) + 1 ^( 33 ) ^( 32 ) + B\ 33) ^'(32)] ^ m ( 3 i) 


+ ^sVu I u m ^(33) S'(32) S m (3i) + ^r*r-'u i u m ^(33) ^(32) ^ m ( 3 i) 


(A.13) 


Using eqs. (A.2) and (A.3) and integrating over u, r and s by eqs. (2.30), (A.5) and 
(A.6), one gets 


lZf(x + ,y + ]y)i =3 = V+U (x + ,y + ; [z(z H 


dz+ / dz£ / dz^ 
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2(k+y 
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A.4 l = 4 term in eq. (2.37) 
Finally, for l = 4, one has 


Kf(x + ,y + -,y)i = 4 = V+ U (x + ,y + - [z(z + )]) f dzf f dz$ 
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(A.15) 


Following the same steps as in the previous sections (and introducing t = u 4 — 113), one 
obtains 

Kf {x + ,y + -,y)i =4 = v + U (x + ,g + ; [z(z + )]) f dz\ [ dz£ f dz£ f dz£ 

J y+ J z+ J z 2 J 
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and, finally, 
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